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Abstract 

/uSR and 7 Li NMR relaxation measurements in frustrated two-dimensional S = 1/2 Heisen- 
berg antiferromagnets on a square lattice are presented. It is found that both in Li2VOSi04 
and Li2VOGe04 , spin dynamics at frequencies well below the Heisenberg exchange frequency are 
present. These dynamics are associated with the motions of walls separating coexisting collinear 
domains with a magnetic wave vector rotated by 90°. 
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Two-dimensional S = 1/2 Heisenberg antiferromagnets (2DQHAF) have attracted a lot 
of interest in the last decade mainly due to their interplay with high T c superconductors^. 
These systems are characterized by strong quantum fluctuations in view of the reduced di- 
mensionality and low spin value, with long range order taking place, in principle, only at 
zero temperature (T). Still it is possible to further enhance quantum fluctuations with- 
out introducing defects, static or itinerant, by frustrating the antiferromagnetic exchange 
interaction. This situation is encountered in the J1-J2 model on a square lattice, where 
the next nearest neighbour interaction J 2 , along the diagonal of the square, competes with 
the nearest neighbour interaction Ji, along the side of the squareBIS!. On increasing J 2 
one expects first a quantum phase transition from a Neel ordered phase to a non-magnetic 
ground state, whose precise nature is still subject of debate!, and then another transition 
to a collinear ground state. The collinear phase can be considered as formed by two in- 
terpenetrating sublattices with a reciprocal orientation of the Neel vectors which classically 
can assume any orientation. However, this degeneracy is lifted by quantum fluctuations 
and, as a consequence of an order from disorder process, just two collinear ground-states are 
realized!: collinear I, with magnetic wave- vector Q = (q x = n/a, q y = 0), and collinear II, 
with Q = (0, ir /a). These two ground states are degenerate and it is not possible a priori to 
tell which of one of them will be realized, unless other interactions are considered. 

Recently two prototypes of frustrated 2DQHAF with J 2 / J\ — 1, Li 2 VOSi04 ( J 2 / J\ — 1.1 
and Ji + J 2 ~ 8.2 K) and Li 2 VOGe0 4 (J2/J1 = 0.9 and J x + J 2 ~ 6 K), have been 
discovered!'!. Li2VOSi04 , which is the most studied of the two compounds, shows a collinear 
ground state always of type I, with the spins pointing along the x direction. The degeneracy 
among the two collinear phases is relieved by a structural distortion occuring just above the 
transition to the collinear phase at T/v!!. For J\ + J 2 > T > T/v the system is characterized 
by a correlated spin dynamics with domains of type I and II extending over a length of the 
order of £, the in-plane correlation length. 

In this manuscript we will provide evidence, based on /zSR and 7 Li NMR measurements, 
that for Ji + J 2 > T > T/v very slow spin dynamics are developed in the two systems, with 
characteristic frequencies u c <C (Ji + J 2 )/cb/^ . These dynamics are associated with the 
motions of domain walls separating the collinear I and II phases which dynamically coexist 
above T/v- New /iSR data showing that both in Li 2 VOSi04 and Li 2 VOGe04 the transition 
to the collinear I ground state is possibly triggered by the in-plane XY anisotropy are also 
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presented. 

The measurements were performed on powder samples prepared by solid state reaction, 
as described in ref. 9. Longitudinal field (LF) and zero field (ZF) /xSR measurements were 
performed at ISIS pulsed source using 29 MeV/c spin-polarized muons. The background 
signal was estimated below T/v, equal to 2.86 K for Li 2 VOSi04 and 2.1 K for Li 2 VOGe04 
, from the amplitude of the slowly decaying oscillation in a transverse field of 100 Gauss. 
For T < T/v, the frequency of the oscillating ZF /zSR signal directly yields the amplitude 
of the local field at the muon site B^, which is proportional to the sublattice magnetization 
(see Fig. 1). The T— dependence of in Li 2 VOGe04 is very close to the one observed in 
Li 2 VOSi0 4 l. In particular, one finds a critical exponent for the sublattice magnetization very 
close to (3 — 0.235, the one predicted for a 2D XY system on a finite siz e§. The absolute value 
of B^, similar for both compounds, allows to estimate where the fi + is probably localized. 
By assuming a simple dipolar coupling with V 4+ magnetic moments, whose absolute value 
should range between 0.6 and 0.25 /iei, two possible /i + sites have been identified. One 
around (Si,Ge)04 tetrahedral oxygens, the other close to VO5 apical oxygen. 

Above T/v, for T < J\ + J 2 , the decay of the muon polarization was observed to be 
rapidly quenched by a LF and to become field-independent above 200 Gauss (see Fig.2). 
The remanent decay, driven by high frequency spin fluctuations, is weakly T-dependent and 
should be compared to NMR spin-lattice relaxation rate l/T^. These observations point 
out that the dominant contribution to the decay of fi + polarization in ZF is not a spin- 
lattice relaxation process induced by field fluctuations at frequencies u> c ~ (Ji + J 2 )/cs//i, 
as we have erroneously put forward on the basis of ZF measurements only!, but it is driven 
by very low-frequency fluctuations with u c <C (J x + J 2 )A;b/^- In this case the decay of 
the polarization is not described by P^{t) oc exp(-Xt), as for LF above 200 Gauss, but 
one has to resort to the complete Kubo-Toyabe functional. In the presence of a LF, if 
the system is not in the purely static or very fast fluctuations regime u c ^> 7V < Ah 2 >, 
with 7 = 2tt x 13.55 kHz/Gauss the muon gyromagnetic ratio and v < Ah 2 > the mean 
square amplitude of field fluctuations, the fit of P^(t) data with the complete Kubo-Toyabe 
function is rather cumbersome00. Nevertheless, if u c > 7V < Ah 2 >, one can resort to the 
analytical expression for P^t) derived by Keren0, which provides a powerful method to fit 
the decay of muon polarization in LF. 

For T ;§> Ji + J 2 a non-negligible T-independent decay, due to nuclear dipolar coupling, 
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is still present. Although this contribution might be disregarded for T — > T/v, in order to 
correctly estimate the frequency of these low-energy dynamics also for T — > Ji + J2, the raw 
/iSR data of have been first divided by Pfj.it) for T ^> J x + J 2 . Then the data have 

been fitted according to 



where the first term describes the spin-lattice relaxation driven by the fast fluctuations, 
while the second one is Keren analytical approximation of Kubo-Toyabe function (see Eqs. 
3 and 4 of Ref. 14), which depends on the intensity of the longitudinal field B, as well as 
on u c and < Ah 2 >. By fixing A (~ 0.015/is -1 ) from the high LF measurements above 
200 Gauss and fitting /xSR data at different magnetic fields, it was possible to estimate 
< Ah 2 > and the T-dependence of r c = 1/lj c (see Fig. 3). We have found that u c is always 
at least 6 times larger than 7V < Ah 2 >, supporting the applicability of Keren analytical 
expression. r c (T) is observed to diverge exponentially on decreasing T, with an effective 
activation barrier around 3 K for Li 2 VOGe04 and close to 8 K for Li 2 VOSi04 . These low 
frequency fluctuations can originate, in principle, either from /i + diffusion or from a slow 
spin dynamics. The first possibility seems quite unlikely if one considers that in perovskites 
the activation energies are usually at least two orders of magnitude largerEl. In particular, in 
Li2VOSiO"4 and Li2VOGe04 the muon is located either close to the apical oxygen of the VO5 
pyramids or close to an oxygen of Si04 tetrahedra where it experiences a strong attractive 
coulomb potential which hinders its diffusion. On the other hand, the second situation is 
more likely once one realizes that the more rapid increase on cooling of r c in Li 2 VOSi04 
seems related to the larger exchange coupling. 

These very slow fluctuations should cause also an enhancement of 7 Li NMR transverse 
relaxation rate I/T2, which was derived from the decay of the echo amplitude after a 7r/2 — 
t — 71 sequence. Two main contributions to I/T2 are present. The first is temperature 
independent and is due to nuclear dipole-dipole interaction among 7 Li nuclei. Its value 
can be estimated from the lattice sum to be (l/T 2 )di P — 5 ms The second one is a 
temperature dependent contribution proportional to the amplitude of the spectral density 
at lu — 0, which sould be very sensitive to these low-frequency fluctuations^. So one can 



P„(t) = exp{-\t)P^t) K {B,u c , < Ah 2 



>) 



(1) 



write that 1/T 2 = (1/T 2 ) dip + (1/T 2 )', with 
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Here h z (t) is the fluctuating component of the local field at 7 Li nuclei along the external 
magnetic field. Now, by assuming oj c ~ (Ji + J 2 )/c£//i, one would derive a contribution to 
1/T 2 of at least one order of magnitude smaller than the experimental one. Moreover, if uj c 
was well above the nuclear Larmor frequency (oul), 1/^2 and 1/Ti should show the same 
T-dependence, while it is evident from Fig. 4 that 1/T 2 starts diverging in a T range where 
1/Ti is still constant. Only a very low-frequency dynamics at u c <C (Ji + ^&b/^ can explain 
the experimental data for 7 Li NMR transverse relaxation rate. Therefore, the temperature 
dependence of (1/T 2 )' oc r c should be characterized by an exponential divergence, with an 
effective activation barrier equal to the one estimated from /xSR relaxation measurements. 
In fact, the barrier derived from 1/T 2 and /iSR decay rate was found to be the same within 
±15 % , pointing out that both techniques are probing the same dynamics and definetely 
ruling out the possibility of fi + diffusion. This implies that there are two main dynamics 
contributing to the spectral density: one with characteristic frequencies of the order of 
(Ji + J2)&b/^, probed by NMR 1/Ti and fiSR decay rate for B > 100 Gauss, the other with 
oj c <C (Ji + J 2 )£;,b//i, which is evident in low field /iSR decay rates and NMR 1/T 2 . 

Since nuclear spin-lattice relaxation rate probes the amplitude of the spectral density 
at lul, one would expect to observe an enhancement in 7 Li 1/Ti on decreasing the Larmor 
frequency towards u c . For this reason 7 Li NMR 1/Ti was measured at B = 3.5 kGauss and 
compared with the one measured at 16 kGauss. One finds that upon cooling, approaching 
T/v, an enhancement of 1/Ti with decreasing field is evident (see Fig. 4), suggesting the 
presence of very low-frequency dynamics. These dynamics are not due to spin diffusion^, 
which in 2D systems gives rise to a logarithmic divergence of 1/Ti. First, because if spin 
diffusion was present one should observe a frequency dependence also for T > 3 K, while 
it is absent (see Fig. 4b). Second, the low-frequency divergence induced by spin-diffusion 
should be cut at frequencies corresponding to energies of the order of the spin anisotropyS, 
which are usually more than three orders of magnitude larger than huj c . 

Now, which is the possible origin of the spin dynamics at frequencies oj c <C (Ji + ^^b/^ ? 
In frustrated 2DQHAF with J\ ~ J 2 the degeneracy of the two collinear ground states leads 
to the coexistence of domains of type I and II above T~ T(T)i, the energy barrier separating 
the two collinear phases. In these vanadates the degeneracy is lifted by a structural distortion 
occurring just above T/v@'i, so that the double potential well describing the energy levels 
becomes asymmetric and the system collapses always in the collinear I phase. On the other 
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hand, for J > T > Tjv, above the structural distortion, the double potential well is symmetric 
and the spectral density is characterized both by a fast dynamic, at oj c ~ ( J\ + ^2)^5/^, and 
a slow thermally activated dynamic with uo c = u exp(—E(T)/T), u being a characteristic 
attempt frequency. These low-frequency fluctuations can originate from the motions of 
domain walls separating the two collinear phases, or, in other terms, from the phase modes 
of soliton strings separating domains^ with a pitch vector Q rotated by 90°. Such a low- 
frequency dynamic is present in one-dimensional QHAF as CuGeOs^, where the slowing 
down of soliton phase modes manifests itself as an increase of the NMR linewidth. uo c can 
be much smaller than (Ji + J 2 )fcB/ft first because of the activated type of dynamics, second 
because 1/uj corresponds to the time required for a domain wall to visit all lattice sites 
within a domain, which increases with the in-plane correlation length £(T) oc exp(27rp s /T), 
p s being an effective spin-stiffnessEl. 

In order to verify if the temperature dependence of r c reported in Fig. 3 for Li 2 VOSi04 
and Li 2 VOGe04 is compatible with such a picture, the data were fitted using for E(T) the 
expression derived by Chandra et al.i 



E(T) 



JlS 2 
2 Jo 



1\ / T 



0.26 - +0.318 



e(T) (3) 



,SJ ' \J 2 S 2 , 

One observes (Fig. 3) that on approaching T N the increase of r c can be well approximated 
by an activated T— dependence, with E(T) given by Eq. 3, with values for 2np s slightly 
below (Ji + J 2 ), as expected for a frustrated systemS. In particular, one finds 2np s ~ 7.6 K 
for Li 2 VOSi04 and 2np s ~ 4.9 K for Li 2 VOGe04 . A departure from this trend is evident for 
T — >• Ji+ J 2 and can be associated with the inadequacy of the simple exponential expression 
used for £(T) at T ~ J\ + J 2 . The consistency of the r c (T) behaviour with the T-dependence 
of E(T) predicted! for the energy barrier separating the two collinear phases, reinforces our 
conclusion in favour of a slow dynamics driven by the motions of domain walls. 

In conclusion, from the analysis of /iSR and 7 Li NMR relaxation rates in Li 2 VOSiO"4 and 
Li 2 VOGe04 it has been shown, for the first time, that in frustrated 2DQHAF, with J\ ~ J 2 , 
low frequency spin dynamics are present at frequencies well below the Heisenberg exchange 
frequency. These dynamics are associated with the motions of domain walls separating the 
collinear I and II phases which coexist above T/v and are characterized by a pitch vector 
rotated by 90°. 
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FIG. 1: Temperature dependence of the local field at the muon site in Li2VOSi04 (circles) and 
Li2VOGeG"4 (squares), normalized by -B^(O) = 313 Gauss and -B^(O) = 339 Gauss, respectively. 
The temperature is normalized by T N = 2.86 K for Li 2 VOSi0 4 and by T/v = 2.1 K for Li 2 VOGe0 4 
. The dotted line shows the critical behaviour expected for a critical exponent (3 = 0.235. 

FIG. 2: a) Time evolution of the muon polarization in Li2VOSi04 in a LF of (from the bottom 
to the top) 5, 10, 20, 50, 900 Gauss, for T = 2.95 K (just above T N ). In b) P^t) for a LF of 5 
Gauss is reported for clarity. The solid lines correspond to the best fits obtained from Eq. 1 with 
7V< A/i 2 > = 0.65 / us- 1 and r c = 0.265 fis. 

FIG. 4: a) Temperature dependence of 7 Li NMR 1/T 2 in Li 2 VOGe0 4 above T N , for H = 3.7 
kGauss. The dotted line shows the best fit according to an activated T— dependence of 1/T 2 with 
an effective barrier of 3.5 K. b) Temperature dependence of 7 Li NMR 1/Ti in Li 2 VOGe04 above 
T/v, for H = 3.7 kGauss (closed circles) and H = 16 kGauss (empty circles). 



FIG. 3: Temperature dependence of the correlation time for the low-frequency fluctuations in 
Li 2 VOSi04 (circles) and in Li 2 VOGe0 4 (squares). These data have been derived from the fit 
of -P M (t) in a LF of 5 gauss according to Eq. 1, after dividing the raw ^SR data by P M (i) for 
T S> J\ + J 2 (see text). The dotted lines show the best fits according to r c = (l/oj )exp(E{T)/T), 
with E(T) given by Eq. 3. The increase in the error bars at high T originates from the increase in 
the weight of the dipolar contribution to P M (t) at higher T. 
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Fig. 3 (P. Carretta et al.) 
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